Kondo Vortices, Zero Modes, and Magnetic Ordering in a Kondo Lattice Model 
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Motivated by the mysteries of the heavy fermion quantum critical point, we investigate the com- 
petition between Kondo screening and magnetic ordering in the honeycomb Kondo lattice at half 
filling. We examine the destruction of the Kondo phase by proliferating vortex configurations in the 
Kondo hybridization order parameter. We find that there are zero modes associated with Kondo 
vortices. Condensing these vortices can lead to the antiferromagnetic phase. 



The quantum critical point (QCP) separating the anti- 
ferromagnetic (AF) phase and the heavy Fermi liquid 
(HFL) phase which is tuned by either magnetic field, 
pressure, or chemical doping is a topic of great current 
interests [1, 2]. The HFL phase is characterized by its 
large Fermi surface, in contrast to the AF phase with a 
small Fermi surface. One suggestion [2. 3] is that the 
Kondo screening and magnetic ordering both collapse at 
the critical point. Clearly this simultaneous collapse is 
the key physics that has to be understood, and it may 
lead to understanding the strong non-Fermi liquid be- 
havior observed in the quantum critical region above the 
QCP. However there are no established theories that ex- 
plain why this simultaneous collapse should happen and 
there are no concrete models in which the simultaneous 
collapse of Kondo and AF scale is studied. 

In this letter we take a step in this direction. Our ap- 
proach is to identify configurations in the hybridization 
mean field theory of the Kondo screening phase which 
may drive the transition to magnetism while at the same 
time kill the hybridization. Natural candidates are vortex 
configurations of the complex mean field hybridization 
amplitude whose structure of vortices may be determined 
within the hybridization mean field theory. Such vortices 
are interesting per se as possible mean field configurations 
in the Kondo screened phase and their core structure 
might offer clues on the transition to magnetism. Sim- 
ilar phenomena happen near deconfined quantum criti- 
cal points in insulating magnets [4]. As a concrete ex- 
ample we study the Kondo-Heisenberg model on a half- 
filled honeycomb lattice. We show that such a vortex 
configuration supports zero modes of the fcrmionic de- 
grees of freedom. This allows the construction of a spin 
triplet vortex. Proliferation of these vortices suppresses 
the Kondo phase and simultaneously produces antiferro- 
magnetic order. Despite this apparently exotic route, in 
this example the critical properties of this transition are 
in the conventional 0(3) universality class. 

The Hamiltonian of the Kondo-Hcisenberg model on 
the honeycomb lattice at half-filling is: 

H^-t^ (cLcja + H.c)j +Jk St-S,+JH S^-Sj, 

{ij}a i (ij) 



where a e {Tii} is the spin flavor. We start inside the 
KI phase. The simplest way to realize the Kondo screen- 
ing phase is to write the localized spin Si at site i by 
employing slave fermions fi'. 



(2) 



where the Hilbert space for the slave fermions is con- 
strained J2afiafia = 1- Thc Kondo Screening phase is 
now realized at the mean field level when the admixture 
bi ~ ifja'^ia) is nonzero and opens a gap in the Dirac 
spectrum [5]. The Heisenberg term also introduces the 
following RVB mixing x^] = {fljjo)- Both 6, - |6,|e'»- 
and Xij = IXij|c*°'^ ^-^e complex- valued. Beyond mean 
field theory, bi is not gauge invariant and never orders 
but the development of the Kondo amphtude \{bi)\ = b 
provides a stiffness which we are going to further ex- 
plore here. In particular we consider proliferating Kondo 
vortex configurations as a way of destroying the Kondo 
coherence. The Kondo field for a KV in the continuum 
limit takes the form b{r) = |6(r)|e^*^^''-' , where \b{r)\ oc r 
as r — ^ and converges to the mean field value b^ as 
r ^ oo. Due to the presence of \{d^ + ia^)b\^ in the 
action, the finite energy configurations are obtained by 
inserting a =F27r gauge flux extended around the vortex 
core. An example of such a conflguration is 



6(r) = feoctanh(r/6)e+'''('-), 
ae(r) = -tanh^(r/A)/r, 



(3) 
(4) 



where the following gauge a{r) = ae{r)6 is chosen for 
the gauge held. We stay with this gauge choice for our 
discussion of the zero mode equations. 

The mean field Kondo Hamiltonian is given by 



H2^-tY + H.c.) 

{ij)a 
{ij)a 



(5) 



where the chemical potentials ^i'^ and /if for the c and 
(1) / fermions have been set to zero. It turns out that the 
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average fermion and electron number is maintained at 
half-filling, even though particle-hole symmetry is de- 
stroyed in the presence of vortices. The proof pro- 
ceeds in two steps. I) The simultaneous transforma- 



tion {bi 



,} converts a vortex to an 



anti-vortex. The wavefunctions become complex conju- 
gated and the density, {c\^Cia) stays the same. II) Un- 
der the additional particle-hole transformation {cm — > 



fly}, the Hamilto- 



nian is invariant. This means that (cJ^Cia) for a vortex 
equals 1 — {claCia) for an antivortex. Combining I and II 
results in (cl^Cia) = 1 — {cl^Cia) = 1/2. The same proof 
goes for 

In the approach we have taken, zero modes play a cru- 
cial role - in their absence, we would have to occupy the 
negative energy states with up and down spins and end 
up with a Dirac sea spin singlet state. To find the zero 
modes we expand the Hamiltonian near the Dirac nodes, 
and set energy=0 to obtain (near the +kB node): 



zT29i -|- iTld2 b 

b* -T2(i9i + ai) - Ti(ic'2 + 02) 



= 0, 



where 7r+ is the column vector [ca+ cb+ fA+ fB+ ) , 
and are Pauli matrices acting on the AB flavors. Oj are 
the cartesian coordinates of a(r) = ae{r)0. The ratio ^ = 
t/x has been eliminated from the zero mode equations 

zero mode equations are studied in the polar coordinates, 
which is natural in our gauge choice. A further rescaling 



after scaling ( '^^^ 



/VC, b VCb. The 



7r+ exp 



transforms this problem to 



ae{p)dp/2 7r+ 



(7) 



d id a{r) 
^ rdd ^ ~ 



CB+ir) + bir)fA+{r) = 0, 



(8) 

which has been studied recently [6], and the existence of 
1 decaying zero mode is shown analytically. The same ar- 
gument goes for the zero mode equations near the —ku 
mode. Therefore counting "f and | spin flavors there are 4 
zero modes in the KV background. We have also modeled 
zero mode differential equations numerically and found 
a decaying zero mode in confirmation of the analytical 
arguments. Furthermore we model the original lattice 
model in a KV background, and in diagonalizing the lat- 
tice Hamiltonian we find the evidence of very small en- 
ergy modes in the Kondo gap which rapidly converge to 
zero as the system size is increased. We have also con- 
firmed that a KV with 27rn flux has 4n zero modes. 
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FIG. 1: 1(a) The schematic figure of the energy levels of the 
Hamiltonian of Eq. (5) in the presence of a KV background. 
An example of KV fields (in the continuum limit) is given in 
Eq. (3), and Eq. (4). The occupied states are colored by blue. 
The Dirac sea in the presence of KV has 2 states less than the 
KI Dirac sea and 2 of the 4 zero modes have to be occupied. 
1(b) An example of a gauge-symmetric (ring) geometry to 
study how v^^ transforms under the 7r/3 rotation around the 
center of the plaquette labeled *. The bold links are the links 
where the gauge field is non-zero. We let the 27r fiux to spread 
over a number of rings A'^a, and choose a symmetric gauge to 
enclose the flux. The numerics is done in the open boundary 
condition. 



After establishing the zero modes, we can now dis- 
cuss the construction of a spin-1 vortex creation oper- 
ator. We define vortex creation operators as operators 
that increase the gauge flux by 27r. If we only limit our- 
selves with states that are connected to MF state, the 27r 
flux-increasing operator contains two terms: 



\aa){bp} 
j-t 



= |G) {DS,-\ 



(G| 



(9) 
(10) 



where a, 6 e {+, — } are nodal, and a, /? £ {t, 1} are spin 
flavors. \G) is the MF ground state, \DS+) is the Dirac 
sea of negative energy states in the presence of the -|-27r 
gauge flux, and is the zero mode annihilation oper- 
ator (with aa flavor) for the state with — 27r gauge flux, 
etc. The Dirac sea in the presence of a KV has 2 states 
less than the KI Dirac sea and therefore 2 of the 4 zero 
modes have to be occupied. There are ways to oc- 
cupy the zero modes. We classify these 6 KV creation 
operators into spin-triplet nodal-singlct, and spin-singlet 
nodal-triplet operators. Since we arc interested in the 
magnetically ordered phases that can arise by condens- 
ing Kondo vortices we focus on the spin 1 and nodal 
singlet operator. The 3 components of this spin triplet 
vortex creation operator is then given by: 



(11) 



where ^ is an arbitrary constant at this point. It is 
straightforward to show that under SU (2) rotation in 
the spin space, v^^ is rotated as a 0(3) vector. 
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TABLE L The table of the transformation of lattice fields un- 
der time-reversal T, charge-conjucation C, and a 7r/3 rotation 
around the center of a plaquette (labeled *) 7?.*y3 . Primed i' 
etc, is just the transformed index under 7?,*y3. Other lattice 
space-group transformations (translations, rotations, and re- 
flections) acts in the same way as in 7?.* ^3 - in that they only 
act on the site indices i — > i'. 



We next wish to see if condensation of the spin triplet 
v'^ can lead to antiferromagnetic order. This requires 

that the v"^ transform identically to the Neel order pa- 
rameter under all symmetries (spin, lattice, time rever- 
sal and charge conjugation). The transformation proper- 
ties of vortex operators under lattice and other symme- 
tries has been addressed recently in the more complicated 
cases of gapless spin liquids. One approach has been to 
measure spin operators expectation values in the vortex 
background after projection [7]. Alternatively the trans- 
formation of each single particle state in the presence of 
a vortex background is studied [8]. The transformation 
of the vortex creation operators can then be read since 
vortex creation operators (in this framework) are con- 
structed from the single particle states. In contrast to the 
projection approach which is a gauge invariant method, 
the second approach only works in gauge-symmetric con- 
figurations. Wc adopt the "single-particle" approach of 
Alicca et at in our problem. In our model we further 
have the advantage that we study the quantum number 
of these vortex excitations in a gapped (KI) phase. 

We first focus on the transformation of under time- 
reversal transformation. Since time-reversal changes +2Tr 
flux to —2tt flux it should send v^^ to v^^ except for pos- 
sible phase factors. Demanding this will require ^ = ±1. 
We sketch the proof here. First we obtain how the in- 
dividual terms in the definition of transform under 



time-reversal: 
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[{icr )a' UfsiifJ. )ab(i(J )aa'{io- )p(3'r)\^a')(bi3') ^^^^ 

= [-(jcr^)(icr2)a*^(ia2)]^^^^ (*M^)afc™^aa')(&/3') 



Similar algebra for the second term in the definition of 



vl^ results in 



(13) 



Comparing this to v^: 

vl = - [a'{ta^)]^^ (il^^)abm^aa)m 



(14) 



i''™(aQ)(f)/3)' 

results in f*^ = 1, therefore ^ = ±1. Due to anti- unitary 
character of the time reversal transformation ± factor 
can not be absorbed by a U(l) rotation. Therefore we 
have two classes of vortex creation operators and their 
condensation lead to very distinct phases. For example 
condensation of can not describe an AF phase as it 
is even under time-reversal. The charge conjugation acts 
similarly to time reversal in that it sends a vortex to an 
antivortex. A careful treatment of charge conjugation 
Tv±. 



results m C : vl_ 



To obtain the transformation of 



under TZ* we 



resort to numerics. We consider a gauge symmetric con- 
figuration [shown in Fig. 1(b)] and apply the method of 
Alicca et al. [8]. Wc find Te^/g : -w^^ The minus 

sign is obtained independent of lattice sizes and vortex 
configurations. This is quite a nontrivial result as all the 
states in the Dirac sea and 2 zero modes contribute to 
this minus sign. 
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TABLE II: The table of the transformation of the 0(3) vortex 
creation operator v^. 



To find how transforms under lattice translations 



02 = 



we use the following identity in the Honeycomb lattice 
7T/3-^a2 = 1' where Oi = (0, V3), and 
(—3/2, —a/3/2) in the units of nearest neighbor 
links. This results in Ta^ ■ v^^ v^^ . The same result 
holds for Ta^ ■ 

Results of the symmetry transformations are summa- 
rized in table II. We see that Re(D_) = + t>l trans- 
forms identically to the standard two sublattice antifer- 
romagnetic Neel order parameter. Thus its condensation 
will lead to the usual Neel order. 

To describe the universality of the resulting magnetic 
phase transition it is convenient to pass to a dual de- 
scription [9, 10] directly in terms of the Kondo vortices. 
As the Kondo hybridization field b is coupled to a gauge 
field, its vortices do not have any long range interactions. 
The dual free energy may then be readily written down 
by demanding invariance under all physical symmetries 
and is given by: 

?=±i (15) 
+s^v1vf + w^\v^ X H y 
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Wc emphasize that, in contrast to the usual boson- vortex 
duaUty, due to the terms here the vorticity is not con- 
served. In other words the free energy is not invariant un- 
der a phase rotation of the vortex fields. This is because 
the gauge field Uij in the original description is compact. 
This allows for instanton configurations where the gauge 
flux can change in units of 27r. However the spin carried 
by the vortices prohibits single instanton events; pairs of 
vortices in a spin singlet can nevertheless be created or 
destroyed as described by the term. 

If r_ < 0, Re(t>_) will condense first, while Im(?;_) 
remains zero and the transition is then described by the 
following free energy 



(16) 



which describes an 0(3) transition for Re(v_). There- 
fore in our theoretical framework, this KV mediated AF 
transition is an 0(3) transition. This is perhaps not un- 
expected, because a charge gap exists on both sides of 
the AF transition and the notion of an onset of Kondo 
screening is an artifact of mean field theory. 
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TABLE III: The transformation of the 8 component fermionic 
field n near the Dirac nodes, fl = IlV^. k', a', and 
are Pauli matrices acting on fermionic 'i'^, spin |J,, nodal -I- — , 
and sublattice AB flavors respectively (see Table IV). 
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TABLE IV: The definition of the 8 component fermionic field 
n. Redundant definitions are ignored. 

The phase characterized by %V(^ x t>| = Re(t>^) x 
Im(?;^) 7^ is also a very interesting phase. This is the 
analogue of the Kane-Mele [11] phase in our model, since 
iv^ X w| transforms the same way as IlfT^K^n (see Ta- 
bles II, III, and IV), and the gap generated by ftfT^K^II 
supports a spin Hall effect. This can be seen by applying 
the external gauge fields A'^^ and A^^ which are coupled 
to charge and spin currents. The effective action for the 
gauge fields is obtained after integrating out the fermions 
SeffiAI^.A"^ = ^ J d^xe^^^xAl^d^Al indicating the spin 
Hall conductivity of magnitude ^ (e = 1) [11, 12]. This 
phase is charactrerized by the spontaneous onset of spin- 
orbit coupling [11]. However to condense Re(t;^) x Im(t>^) 
we should tune both t^ and rj, and the transition from KI 



to this Kane-Mele phase is controlled by a multi-critical 
point. 

The KV excitation wc considered is a charge-0 spin-1 
excitation and deep in the KI phase is not distinguish- 
able from a particle- hole excitation which is accompanied 
by a spin-fiip. However wc have checked the energetics 
of the KV excitation and have found that the KV ex- 
citation costs less energy than a simple KI particle-hole 
excitation. Therefore in the band theory terminology (as 
far as long distance effects of these excitations are con- 
cerned) the KV (in our model) may correspond to triplet 
excitonic bound states in the KI gap and the onset of 
AF may be viewed as the condensation of these triplet 
excitons. 

In summary we go beyond mean field theory to de- 
scribe the destruction of the Kondo admixture parame- 
ter due to the proliferation of Kondo vertices. Wc find 
that the vortex field carries spin 1 quantum numbers and 
its condensation leads to AF ordering. In this way the 
suppression of Kondo screening is linked to the appear- 
ance of AF. The drawback of our model is that a charge 
gap always exist, both in the KI and in the AF phase 
and is nonzero at the onset of AF. As a result there is no 
sharp onset for the KI, in contrast to the HFL case where 
the Knodo phase is characterized by the sharp onset of 
a large Fermi surface. Nevertheless, we hope that this 
line of investigation might shed some light on the more 
difficult problem of the AF-HFL quantum critical point. 
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